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The chiral limit of finite- volume QCD is the e-regime of the theory. We discuss how this regime can be used for 
determining low-energy observables of QCD by means of comparisons between lattice simulations and quenched 
and unquenched chiral perturbation theory. The quenched theory suffers in the e-regime from "quenched finite 
volume logs", the finite- volume analogs of quenched chiral logs. 



When the chiral limit m — ► of QCD is taken 
at finite four-volume V , one will eventually, no 
matter how big V is, reach a regime where the 
pion Compton wavelength l/m n exceeds the lin- 
ear size L = V 1 / 4 of the box. If the volume at the 
same time is much bigger than the QCD scale, one 
is in an interesting region, the e-regime of QCD, 
where the spontaneously broken chiral symmetry 
is dictating the interactions of the pions, while, 
simultaneously, these pions do not even fit in- 
side the box |1I2| . It is important that there is 
a huge separation of scales here: Seen from the 
mass scale of the pseudo-Goldstone bosons m,, 
the rest of hadronic physics is effectively occur- 
ing near the cut-off. This isolates a neat little 
corner of QCD which can be studied by analyt- 
ical means, and the results are completely non- 
perturbative in the strong coupling constant. Far 
from being an academic playground, this e-regime 
can thus be very useful for lattice gauge theory. It 
is a new regime in which to do numerical compu- 
tations in lattice QCD, and where it is extremely 
important to get a precise handle on the lowest- 
lying Dirac eigenvalues PJ- These smallest eigen- 
values dominate physical low-energy observables 
in this regime because the mass is taken to zero. 

The term e-expansion is used to denote a re- 
ordered chiral perturbation theory. Normally, 
for volumes infinite or very large compared with 
the pion Compton wavelength, the expansion pa- 
rameter is basically "momentum", i.e. p 2 ~ 
(hence the alternative name "p-expansion" ) . In 
the e-regime such an expansion is doomed to fail 
due to the enhanced role played by momentum 
zero-modes. For these zero-modes the pion prop- 
agator would blow up, were it not for the mass 
term mj. As the chiral limit is taken at fixed four- 
volume, even this term will blow up. Because the 
problem only concerns zero-modes of momentum, 
a cure is fortunately available The idea is to 



do the path integral of the zero-momentum modes 
exactly, and treating the remaining modes per- 
turbatively. This is possible, since the non-zero 
modes will have a gap on the order of l/L down 
to the zero modes. The small quantity e in which 
one is expanding is precisely l/L. There is a fasci- 
nating relation between the pertinent zero-mode 
integral of the chiral Lagrangian and chiral Ran- 
dom Matrix Theory. Much of the recent progress 
in understanding the e-regime of QCD stems from 
this relation to universal Random Matrix Theory 
results |4I5| . 

While this e-expansion works well for full QCD 
in the chiral limit jH], fully quenched analogs of 
chiral perturbation theory |7I8| suffer also in the 
e-regime from the double-pole effects of the sin- 
glet field |9lll)j . The probably most drastic effect 
of quenching occurs already at the one-loop cor- 
rection to the quenched chiral condensate. Just 
as quenched chiral perturbation theory in the p- 
expansion shows internal inconsistency by gen- 
erating arbitrarily large corrections in the form 
of quenched chiral logs [711 1| . quenched chiral 
perturbation theory in the e-expansion produces 
"quenched finite volume logs" [2}. Although the 
e-expansion should become better and better as 
the volume V is increased, a logarithmic correc- 
tion prevents this limit from being taken! The 
origin of this trouble lies in the fact that the 
quenched e-expansion has one additional expan- 
sion parameter, mo, the mass term of the flavor 
singlet. Contrary to l/L, this parameter is fixed 
by physics, and cannot be tuned. It means that 
the correction term formally can become larger 
than the leading term, indicating that the whole 
expansion, at least as done in this straightforward 
way, is sick. The extent to which one can ignore 
this problem of quenched QCD will be discussed 
below. Fortunately full and, to some extent, par- 
tially quenched QCD does not have this diffi- 
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culty, and there the e-expansion converges nicely 
as V — > oo. 

For simplicity partially or fully quenched chi- 
ral perturbation theory will be discussed in the 
framework of the replica formalism |8ll(Jj . All 
results reviewed here have actually been calcu- 
lated in the supersymmetric formulation [7] as 
well, as a check. With replicas, the full or par- 
tially quenched limit in chiral perturbation the- 
ory is quite simple, and differs only in the fully 
quenched case from ordinary chiral perturbation 
theory, and then only in the necessity of includ- 
ing the flavor-singlet field together with the usual 
Goldstone bosons. The flavor singlet cannot be 
removed in the fully quenched theory because 
there is no replica limit N — > of the flavor group 
SU(iV). As in the supersymmetric formulation, 
one "regularizes" by adding a massive singlet field 
<&q(x), whose mass mo in the replica formulation 
serves as a deformation of the symmetry group 
away from SU(A). The SU(iV) limit is reached 
as mo — > oo, but that is a limit that cannot be 
taken before N — > 0. The chiral Lagrangian then 
reads 
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Here the ellipses denote terms that are sub- 
leading in both the p and e expansions. The 
a-term in Q is often omitted (because it has a 
power of p 2 compared to the mo-term), but can 
in any case easily be reinstated in all formulas in 
momentum space by m§ — > m§ + ap 2 . The con- 
stant £ is the (infinite volume) chiral condensate, 
and is the pion decay constant. The vacuum 
angle 6 is introduced only in order to be able to 
project on sectors of fixed topological charge v 
(see the second of ref. £Q). 

When ordinary chiral perturbation theory 
based on breaks down in the chiral limit at 
finite volume it happens because the zero mo- 
mentum modes cannot be treated perturbatively. 
The cure is to treat just these modes exactly. In 



practice this is done by factoring 

U(x) = U Q exv[V2iZ{x)/F^] (2) 

where Uo is the zero-mode part, and all flutuation 
fields are included in £(x). In conventional chiral 
perturbation theory it suffices to expand also Uq 
around the origin; here it must be allowed to fluc- 
tuate all over the Goldstone manifold, and the in- 
tegral over Uo will be evaluated exactly. The full 
integral to be performed is not just over the ac- 
tion JJJ, but with the observables inserted in the 
path integral as well. In practice this is done by 
evaluating a suitable generating function of the 
zero modes exactly. 

The quenched finite-volume logarithm 

Consider the (mass-dependent) quenched chiral 
condensate 



£„(m) = lira -^— r ^—\nZ u 

n^q NV dm 



(3) 



where Z v is the partition function in a sector of 
topological charge v. To leading order in the e- 
expansion this can be evaluated analytically on 
the basis of the zero momentum integral over Z7o 
0, 
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where /i = mY>V , and / and K are modified 
Bcsscl functions. This leading order expression 
has been compared to lattice gauge theory data 
with quite nice agreement |12j . But an unpleas- 
ant result emerges when we evaluate the one- 
loop correction to this result. In the quenched 
e-expansion this can be done by computing the 
one-loop ^-integral saturation of 



{U + Ul) dx £(x) 



1 / mX 
hm — ( „„ Tr 

n^o N \ 2F 2 



with, in the replica formalism, N replica quarks 
of equal mass m (the quenched pseudo-Goldstone 
bosons thus have common mass M = 2mH/F 2 ). 
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To evaluate it, we need the propagators of, in 
the quark basis, the off-diagonal mesons $y ~ 



P 2 +M 2 ' 



(5) 



and the more complicated propagator for the di- 
agonal combination $>n ~ tpitpi, 



p 2 + M 2 {p 2 + M 2 ) 2 T{p 2 ) ' 



where |S] 



2s- 1 , (mg+ap 2 )AT 
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This correction can be rewritten in precisely the 
same form as the original mass term in the La- 
grangian. We can thus re-exponentiate, and find 
that to this order in the quenched e-expansion we 
can use the same effective Lagrangian, but with a 
shifted mass term. In detail, we replace \x = rnSV 
by fi', where 
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(11) 



(and hence equals unity in the quenched N — > 
limit). Because the integration over £ excludes 
zero modes, it is convenient to define 



e ipx e ipx 



(8) 



and note that these are just obtained from the 
second and third terms in the Taylor expansion 
of the pion propagator at x = 0: 



V ^ 



a ipx 



V ^ p 2 + M 2 M 2 V 

v 



+ AO) - M 2 G{x) , (9) 



in an obvious notation. Performing the contrac- 
tions, one gets, successively, 



and compute the one-loop improved chiral con- 
densate from 



J? — lim — — — - InZ(a') ■ — 
E jv-o N dfj,' v ' n 



(12) 



Alternatively, to this order one is effectively work- 
ing with a shifted ^-parameter in finite volume. 
To bring out the volume dependence it pays to 
rewrite it as 



1 



E N r F 2 



= 1 + F^KG(0) + aA(0)], (13) 



valid to this order and in the chiral limit. In di- 
mensional rcgularization A(0) is finite |13j. 



A(0) 



(14) 



1 VmT, , u 



(iV-l)A(M 2 ) + i^G(M 2 )l 



1 VmY, , 4. 
= lim — — r Tr {U + W x 
w^o AT 2F 2 

{(JV- 1)A(M 2 ) + A(M 2 ) 



where j3\ is a universal "shape coefficient" that 
depends only on the geometry of the finite vol- 
ume. There is, even in dimensional regular- 
ization, an ultaviolet divergence associated with 
G(0), 



G(0) = & + JL(ln(LAo)-ci) , 



(15) 
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where /3 2 is another universal shape coefficient 
|13| . A is the momentum space subtraction 
point, and c\ is, still in dimensional regulariza- 
tion, 



d = -J__i(r'(l) + H-ln(47r) + 0(d-4)) .(16) 



In fact, this divergence found here in the e- 
expansion matches precisely the new quenched 
one-loop counterterm of the infinite-volume the- 
ory found by Colangelo and Pallante This 
is as it should be: the divergence structure at fi- 
nite volume is just as in infinite volume, no new 
counterterms are needed, and the already exist- 
ing infinite- volume counterterms are cancelled by 
similar divergences in the finite-volume theory. 

If we plug HI4fl and l|I5|l into the expression for 
the effective chiral condensate parameter S e // in 
(|13|) we find a logarithmic divergence in L, the 
quenched finite volume logarithm. As this is sup- 
posed to be an expansion in 1/L, something has 
obviously gone terribly wrong. It is worth not- 
ing that the logarithmic term (|15|) occurs even in 
the unqucnchcd theory, but there always accom- 
panied by inverse powers of L, which makes its 
appearance there innocuous. We trace the origin 
of the logarithmic disaster in the quenched theory 
back to the strong infrared divergence associated 
with a double-pole propagator. This failure of the 
quenched e-expansion occurs because there really 
is one more expansion parameter involved, mo, 
which is not a tunable quantity. Rather mo is re- 
lated to the flavor singlet mass of the unquenched 
theory, and it remains a finite number even in 
the chiral limit. The quenched finite- volume log- 
arithm occurs in the two other chiral symmetry 
breaking classes too ^HJ) and thus appears to be 
unavoidable consequence of quenching. The fact 
that what should be a small correction to the 
leading-order result can become big means that 
the expansion as it stands cannot be trusted. To 
this order, the best one can do is to view it as a 
prediction for the finite-volume scaling behavior 
of nearby volumes V\ and V<j. Then the pertur- 
bation theory to this order still makes sense, and 



one gets the prediction 
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to this order. Taken at face value it suggests 
that the infinite volume chiral condensate one at- 
tempts to extract from the quenched theory may 
keep growing as the volume is increased. This 
type of behavior has often been suspected, and 
there is also numerical evidence from two space- 
time dimensions to support it |I6) . 

Quenched Correlation Functions 

The e-expansion gives rise to a new chiral per- 
turbation theory for the low-energy observables, 
and in particular n-point functions of mesonic op- 
erators can be worked out in this framework. The 
quenched finite volume logarithms, and quench- 
ing artifacts in general, will put restrictions on the 
applicability of these expressions in the quenched 
theory. Nevertheless, since this is where the main 
numerical focus is with present-day techniques for 
fermion simulations with good chiral properties, 
it is of interest to determine these quenched ex- 
pressions. So far computations have been per- 
formed for two-point functions of both scalars 
(S) and psuedoscalars (P) JU| as well as vec- 
tors (V) and axial vectors (A) J7j. Also certain 
quenched three-point functions have been com- 
puted JSj- Until this year only one preliminary 
numerical study using domain wall fermions had 
explored these quenched correlation functions in 
the e- regime JH] , but one month ago new overlap- 
fcrmion results were reported |2Q| . 

For the flavor singlet correlation functions there 
is no need to work with anything but one valence 
quark, N v = 1. The relevant operators at quark 
level read, in a replica notation, 



S°(x) 
P°(x) 



(18) 



where the N x N matrix is unity in the first 
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entry, and zero elsewhere. The quenched limit is 
taken by sending TV — > ^U]- Sources s(x) and 
p(x) for such operators are readily transcribed 
into the effective Lagrangian by means of replac- 
ing the source (mass matrix) 



M — *• M + s(x)I Nv + ip(x)I N „ 



(19) 



in the chiral Lagrangian 0J. Flavored sources 
are added analogously: we now take N v to be 
(at least) 2, define Jjv„ correspondingly, and then 
transcribe the quark currents 



S a (x) EE <lp(x)t a I Nv 1p(x) 
P a {x) EE lPit a j 5 I N ^{ X ) 



(20) 



into the effective Lagrangian by means of the 
sources 



M M + s a (x)t a I Nv + ip a {x)t a I N „ . 



(21) 



Similarly for the vector and axial vector currents 
which at the quark level read 



TO 
a;(x) 



ip(x)i'y^t a I Nv ip(x) 

4>{x)ijn75t a lNv4>(x) 



(22) 



and which can be assigned sources in the effective 
theory through a replacement of d i JJ{x) by a co- 
variant derivative, as in ordinary chiral perturba- 
tion theory. In the e-expansion the pertinent cor- 
relation functions are then evaluated by perform- 
ing the loop epxansion in the fluctuating £(x), 
while performing the zero-mode integral over Uq 
exactly. It is this last step which is the tricky 
part. Fortunately all analytical results required 
for N v = 2 compuations are available in the liter- 
ature after suitable manipulations |10|. 

The analytical results for these two-point func- 
tions are fairly simple in structure, and all depend 
on two kinematical functions 



hi(r) 



d 3 x A(x) 



h 2 {T) 



1 

24 



K 2' 12 

d 3 x G(x) 



(23) 



where T is the extent in the "time" direction and 
t ee t/T. These two functions replace the simple 
cosh-function of the usual two-point functions in 
conventional chiral perturbation theory. 

To give an example of the kind of expression 
one gets, consider the sum ^Oj 

(S a (x)S a {0)) + (P a (x)P a (0)) = A(ji)+ 

|- [Afc) - B(p)(mlG(x) + aA(x))] , (24) 

where A(/i) and B([i) are known explicit func- 
tions of /j,: 



m 

E 
B(ii) 



A* 



j_ ( Kig l E„(/i) 

N c \ E n E 



(25) 



The function Ej,(/i) is as given in eq. |@}, and the 
"1-loop improvement" is as dictated by the 1-loop 
correction to the chiral condensate as discussed 
above, 



K- loop (v) = W)- ■ 



(26) 



Taking, as usual on the lattice, the projection 
onto zero momentum of this sum of correlation 
functions is then straightforward and expresses 
the spatial integral of this sum in terms of the 
kinematical functions hi(r) and h<z{T). Although 
this looks quite different from conventional two- 
point functions on the lattice, one can check ex- 
plicitly that the above result matches the expres- 
sion from ordinary chiral perturbation theory in 
the region of overlap. The relation ^ is crucial 
in establishing this. 
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We have already mentioned that in the e-regime 
there is strong dependence on gauge held topol- 
ogy. This is clearly seen in the correlation func- 
tions which vary strongly with increasing v. Ex- 
ceptions are suitable combinations of correlation 
functions such as the one discussed above, where 
the quenched poles at zero quark mass in non- 
trivial gauge field backgrounds cancel at any fixed 
v 7^ 0. Zero topology, v = 0, is yet again special: 
there is no quark mass pole but still a strong in- 
frared sensitivity to the quark mass, see ref. ^U] 
for details. 

Even when explicit power-like divergences in 
the inverse quark mass are absent, these quenched 
correlation functions all suffer from the quenched 
finite volume logs. This is clear from the example 
(S a {x)S a {0) + P a {x)P a {0)) discussed above since 
the 1-loop improved chiral condensate E 1_ioop (/i) 
enters directly in this correlation function. Again, 
taken at face value the analytical prediction is a 
divergence in the correlation function. But the 
divergence is also here not to be trusted, since 
it arises as a correction to the leading order re- 
sult. As for the chiral condensate itself one can 
hope that the analytical prediction gives the scal- 
ing between two nearby volumes, where no diver- 
gences arise. It should also be stressed that these 
quenched finite logs are by now means worse than 
the more familiar quenched chiral logs; only the 
fact that in the e-regime one really aims at taking 
the chiral limit makes the issue of these quenched 
infrared divergences more urgent to resolve. 

If one computes vector and axial vector corre- 
lation functions, the quenched finite volume logs 
do not appear neither at leading nor subleading 
order in the e-expansion |17|. This should not 
be taken as evidence that these correlation func- 
tions are free of these logs, only that they do 
not appear at these first 2 non-trivial orders (cor- 
rections to the quenched chiral condensate, and 
hence quenched finite logs are almost bound to 
appear eventually). Actually one of these, the 
quenched vector- vector correlation is free of these 
logs, but in a trivial way: 



to all orders in the quenched e-expansion |17| ! 
The proof of this statement can be found in ref. 
|17j (it is a peculiar quenched artifact with a sim- 
ple diagrammatic interpretation). The 2-point 
function of axial vectors also shows a quite re- 
markable behavior. Up to and including next-to- 
leading order, 

J d 3 x{A « (x)A» (0))=^ (l+ 4mS ;^ )r2 Mrj) (28) 

This result comes about after a series of non- 
trivial cancellations where all double-pole terms 
in the quenched e-expansion have cancelled each 
other. To leading order this gives a very simple 
prediction that does not refer to quenched pa- 
rameters at all: the axial vector 2-point function 
directly measures the pion decay constant F w . In- 
cluding next-to-leading order it also determines 
the quenched chiral condensate (but actual nu- 
merical simulations may have difficulty disentan- 
gling this correction from the leading result |2l)|). 

Predictions for full QCD: 

The difficulties of the quenched e-expansion 
are fortunately absent in the full theory. The 
quenched finite volums logs become, in the full 
theory, multiplied by inverse powers of the vol- 
ume, and are hence rendered harmless. This is 
entirely analogous to the quenched chiral logs 
becoming replaced by ordinary (harmless) chi- 
ral logs in full QCD in the p-expansion. Just as 
ordinary chiral logs are multiplied by powers of 
the pion mass, the finite volume logs in the e- 
expansion become multiplied by inverse powers 
of the volume. 

Analytical predictions for full QCD in the e- 
regime was worked out long ago by Hansen 
and they have recently been re-done |10I17| for 
sectors of fixed gauge field topology. These pre- 
dictions in sector of fixed topology may eventually 
may be of most interest from a numerical point 
of view, in particular since there are distinct pre- 
dictions for each particular ^-sector. 



(V a (x)V a (0)) = 



(27) 



Also in the full theory it turns out that all 
zero-mode integrals can be expressed in terms of 
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the mass-dependent chiral condensate £j,(/Lt) and 
derivatives thereof. As in the quenched theory it 
pays to define, to one-loop order a shifted mass 
which here turns out to be £Q 



Ml 
/' 



N f 



1 Pi 



N f F 2 VV 



(29) 



All analytical predictions are then expressed in 
terms of the kinematical functions hi(r) and 
hstr) up to an including next-to-leading order in 
the e-expansion (additional numerical factors re- 
lated to the lattice volume also appear). Details 
can be found in refs |1(JI17| . 
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